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Abstract 

We have analyzed 116 fb~^ of data collected by the BABAR detector for S"^ production. In this 
paper we describe the observation of production from cc continuum and from B decays, with 
the S^} decaying into S~tt'^ and 17~X+ modes. The ratio of the branching fractions of the 
decays into these two hnal states measured in the cc continuum is: 

Rf "0 _^ []-K+) 

— -f = 0.296 ± 0.018 (stat.) ± 0.030 (sys.). 
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1 INTRODUCTION 



Little is known about charmed baryons today even though decades have passed since the discovery 
of charm. The high-luminosity 5-factories present excellent opportunities to study the production 
and decay of charmed baryons with high precision. 

We present the observation of the (csd)^ charmed baryon in two decay modes: 

"c ^ " vr . 

The ratio of branching fractions of these decay modes has been predicted to be approximately 
0.32 ^ using a spectator quark model calculation. This figure has a substantial theoretical un- 
certainty. It has been measured previously by the CLEO collaboration; their result was consistent 
with this prediction but had a large statistical uncertainty |2]. 

We measure the ratio of the branching fractions of — > Q~K^ and S~'k~^ from the 

continuum production of e^e~ — > cc, where the hyperons are reconstructed through the following 
decays: 

Q- ^ AK~ 

A ^ P TT~ . 

Since the two final states are topologically similar, quite a few systematic uncertainties cancel in 
the ratio of the branching fractions. We also observe signals for T(4S') BB -|- X in both 

final states, where X represents the rest of the event. Although copious production of S"^ and S'^ 
in B decays has been predicted [Hj, this process has been observed previously only by CLEO, with 
a significance of approximately 3 cj in the — s- S^-k^ decay mode and approximately 4 cr in a 
related S'^ decay mode 0]. 

2 THE BABm DETECTOR AND DATASET 

The data for this analysis are collected with the BABAR detector at the PEP-II asymmetric e^e~ 
collider; a total integrated luminosity of 116 fb~^ is used. A five-layer silicon vertex detector (SVT) 
and a 40-layer drift chamber (DCH) form the tracking system. The drift chamber is surrounded 
by the DIRC, a detector of internally reflected Cherenkov light, which provides additional charged 
particle identification (PID). These are enclosed in a CsI(Tl) electromagnetic calorimeter (EMC). 
The detector assembly is embedded in a 1.5 T superconducting magnet. Further details of the 
BABAR detector are given elsewhere [H]. 

The data collected are processed through the standard BABAR reconstruction software. In the 
present analysis, we use an integrated luminosity of 105.4 fb^^ collected at the T{AS) resonance, 
and 10.7 fb~^ collected below the T(45) threshold. We refer to these as the on-peak and the 
off-peak data samples, respectively. 

®A11 channels imply the charge conjugates as well, unless otherwise specified. 
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Figure 1: Invariant masses of reconstructed and preselected (a) S and (b) Q candidates from 
subsamples of data. 



3 ANALYSIS METHOD 

The two decay modes Q~ — > AK~ and S~ Att~ are topologically similar. The hyperons involved 
are long-lived. The decays close to the production vertex: cr = 34^*^2 A*™ 0- 



3.1 Selection of Events 

The S^! reconstruction takes place in three stages: 

• pre-selection of events containing a A, 

• pre-selection of events containing either S~ or i7~ from the A sample, and 

• construction of candidates. 

The A is reconstructed by identifying a proton and combining it with an oppositely charged 
track. A candidates within a 3a {a is the fitted mass resolution) range of the central value are 
then used for reconstruction of S~ and J7~ by vertexing it with a negatively charged track; and 
the A mass is constrained at the nominal value [Sj. For f2^ reconstruction, the K~ is required to 
be identified as a kaon. 

To improve signal-to-noise ratio, a minimum decay distance of 2.5 mm between the primary 
vertex and the S~ decay vertex in the plane perpendicular to the beam direction is required; for 
the f2~ , the required distance is 1.5 mm. In addition, the "signed" flight distance^ between the A 
and the f2~ decay vertex is required to be at least 3 mm. 

Figures ^ (a) and (b) show invariant mass distributions of Ait~ and AK~ respectively, from 
subsamples of the data. Superimposed on the plots are fits to a double Gaussian (single Gaussian) 
for the together with a linear background. The fitted masses and resolutions of data and 

Monte Carlo are consistent within known systematic effects. 

"signed" flight length is where the displacement and the momentum vector of the particle are required to be 
less than 90° apart, i.e., p ■ r > 0, where r denotes the distance from the production point to the decay point of 
particle X in the xy-plane. 
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Table 1: Fit results for 



Yield in 


S-7T+ 


[2-K+ 


On-peak Data 


7614 ± 545 


906 ± 54 


Off-peak Data 


450 ± 39 


78 ± 10 


On- and off-peak 
Data, p* > 1.8 GeV/c, 
in COS0* range: 


4058 ± 319 
(-0.8 < cos 6** < 0.8) 


655 ± 43 
(-0.8 < COS0* < 0.6) 



3.2 Reconstruction 

The selection criteria are finalized using subsamples of the data as a precaution against selection 
bias. The subsamples used are of size 20 fb-i and 40 fb"^ for the S-tt+ and n-K+ modes 
respectively. The final results are obtained using the entire 116 fb~^ sample, including these 
subsamples. In each case a Scr mass range around the central value is used. 

Each resulting S~ candidate is then vertexed with an oppositely charged pion for the S^tt^ final 
state. Likewise, each f}~ candidate is vertexed with a positively charged track identified as a kaon 
for the f2~K~^ final state. The resulting invariant mass distributions for the candidates from 
the on-peak data sample are shown in Figures |2l (a) and (b) for S^tt^ and Q^K^ combinations, 
respectively. The mass distributions from the off-peak data sample are shown in Figures [21 (c) 
and (d) again for E~'k^ and f2~ combinations, respectively. A clear peak is evident in all 
four spectra. The fitted distributions are superimposed on the plots. In each case we use a single 
Gaussian shape on a linear background, except for (b) in which a much better fit is obtained by 
using a double Gaussian shape on a linear background. The fit results are listed in Tabled 

3.3 Simulation 

Events corresponding to the e+e~ — > cc — > S")? -|- X process are generated, with the decays into 
the two desired decay modes. PYTHIA 7 is used for the cc fragmentation and GEANT4 |H] is used 
to simulate the detector response. These events are then reconstructed and the selection criteria 
applied. Samples of 90,000 events for the S~ti~^ final state and 60,000 events for the Q~K^ final 
state are generated. To investigate possible background contributions, generic e"'"e~ — > qq {u, d, s, c} 
continuum Monte Carlo events are processed through the complete analysis program sequence. 
The e~^e~ — > cc sample corresponds to an integrated luminosity of 64 fb~^, and the combined 
uu,dd,ss sample corresponds to 33 fb^^. In addition, 22,000 events are generated according to 
T(4S') BB + X, and processed through the complete analysis chain. 

3.4 Background Contributions 

We analyze the generic cc Monte Carlo events. No evidence of a peaking background is observed. 
The distribution of the background events can be fitted with a linear shape in the S^tt^ channel. 
For the Q~K^ channel, the distribution of the reconstructed events is flat. The events from generic 
uu, dd, s's Monte Carlo do not show any peaking either. 
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Figure 2: The invariant mass distributions for candidates, shown for (a) ^ vr"*" in on-peak data, 
(b) Q~K~^ in on-peak data, (c) H~7r+ in off-peak data, and (d) Q~K^ in off-peak data. 
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P* (GeV/c) 



Figure 3: p* distributions from reconstructed in Monte Carlo. The dashed red line shows 
produced in B decays, and the solid blue line shows cc production of S^. Normalizations are 
arbitrary. No background is present. 

3.5 5"° Production from cc Continuum and from T{AS) BB ^ + X 

The B ^ S^ + X Monte Carlo events are instructive in separating the contribution originating 
from B decays from those originating from the cc continuum production. Figure |21 shows the 
distribution of the momentum of the reconstructed S'^'s in the center-of-mass frame {p*) from the 
Monte Carlo®. These are also "truth-matched", i.e., where the reconstructed information matches 
the generated information. 

The dashed red line shows the p* distribution of S'^'s originating from B decays and the solid 
blue line shows that from cc continuum. The normalizations in this figure are arbitrary. The p* 
distribution from B decays does not extend beyond 2 GeV/c, purely from kinematics, whereas the 
distribution from the continuum peaks at much higher p* values. 

The p* distributions from the signal regions for the off-peak data are shown without any 
efficiency correction in Figures0](a) and (b), for the S~ti^ and the Q~ final states, respectively. 
These data are collected below hb production threshold, and therefore represent production from 
continuum only. The background under the signal in the data is estimated from the sidebands 
in the reconstructed mass spectrum and then removed. These p* distributions, peaked around 
3 GeV/c, clearly indicate production from cc continuum. 

Figures [S] (a) and (b) show the p* distribution in on-peak data from the candidates in 
the S~TT~^ and Q~K^ final states, respectively, after background subtraction, again without any 
efficiency correction. The peaks below 1.5 GeV/c in both plots clearly represent production 
from B decays, as evident from Figures |31 and 0] 

3.6 Analysis and Efficiency Correction for cc — > ° + X 

For further analysis the on- and off-peak data samples are combined; to isolate the cc production of 
events with p* > l.S GeV/c are selected. In order to avoid large fluctuations from the edges of 
the phase-space and detector acceptance effects, we also require —0.8 < cos 6** < 0.8 for the S~tt^ 

*The decay of the B into S'^ is modelled using PYTHIA 13 fragmentation. 
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Figure 4: Sideband-subtracted distribution of reconstructed candidates in off-peak data 
without efficiency correction in (a) S~tt'^ and (b) Q^K^ mode. Most of the signal is produced at 
higher p* as expected. 
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Figure 5: Sideband-subtracted p* distribution from candidates in on-peak data without effi- 
ciency correction in (a) S'^vr^ and (b) Q^K^ mode. The lower peak below p* < 1.5 GeV/c is 
primarily from the production from B decays as evident from Figure |31 and Figure [l] 
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m(S-7t*) (GeV/c') m(n-K*) (GeV/c^) 

(a) (b) 

Figure 6: Invariant mass spectra for with p* > 1.8 GeV/c from (a) S'^vr^ final state with 
-0.8 < cose* < 0.8, and (b) Q-R^ final state with -0.8 < cos 6** < 0.6 using on- and off-peak 
data. Both spectra are fitted with a double Gaussian for signal and a linear background shape. 



mode and —0.8 < cos 9* < 0.6 the for f2~K~^ mode, where 6* is the polar angle of the candidate 
with respect to the collision axis in the center-of-mass frame. 

Figures El (a) and (b) show the invariant mass spectra in these cos 9* ranges with p* > 
1.8 GeV/c for the combined on- and off-peak data sample, for S~Tr^ and Q^K^ decay modes, 
respectively. The fit results are presented in Tabled 

The efficiency is calculated from signal Monte Carlo events as a function of p* and cos0* in 
each of the two decay modes. For each decay mode, a fifteen-parameter two-dimensional fit gives 
a smooth parameterization of the efficiency with small statistical unceratinty. We then correct the 
data distribution by weighting each event in the spectrum inversely by its efficiency according to 
the event's position in (p*, cos^*) space. After correcting for efficiency we obtain 19375±393 events 
in the E~it~^ mode and 4866 ± 283 events in the fi~ mode^. 

The angular distribution of the data is well-described by a (1 -|- cos^ 9*) function; we use this to 
estimate the fractions of the which are expected to lie in the selected angular regions for each 
mode. Extrapolating from the fiducial region into the full range —1 < cos9* < +1, we obtain the 
total numbers of signal events for p > 1.8 GeV/c are 26621 ±540 and 7874 it 458 for the S~Tr~^ and 
Q^'K^ modes, respectively. We thus obtain the ratio of branching fractions: 

—77^ „_ ' = 0.296 ± 0.018 (stat.). 

4 STUDY OF SYSTEMATIC UNCERTAINTIES 

We evaluate several sources of systematic uncertainties, described below and summarized in TablelS 
Adding all of these uncertainties^'' in quadrature, we obtain a total absolute systematic uncertainty 
of 0.030 on the ratio of the branching fractions. 



®The A branching fraction is not taken into account in these numbers, since this cancels in the ratio 
"No be 
neglected. 



^"No baryon polarization is considered in the present analysis and any systematic uncertainty due to this is 
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• We vary the signal shape, the background shape, and the fit range, and use a simple event- 
counting method. From the deviations observed, we assign systematic uncertainties for the 
use of a binned fit and for the particular technique used, adding them in quadrature. 

• We repeat the analysis with (a) a parameterization of the efficiency with a similar function 
with nine parameters, and (b) a simple efficiency calculation in two dimensional bins from 
Monte Carlo. The discrepancy observed between the main result and the result from (a) is 
assigned as the systematic uncertainty. 

• We vary the range in cos 6* used for the two final states. We also vary the cos 9* distribution 
used for the extrapolation. The combined systematic uncertainty is taken to be the sum in 
quadrature of the uncertainty due to using different cos 9* ranges for the two modes and the 
uncertainty due to the choice of extrapolation function. 

• We take into account an uncertainty due to the finite size of the Monte Carlo sample used to 
estimate the efficiencies. 

• Approximately 1% of selected events contain multiple candidates in the signal range 
with one or more tracks in common. We retain all such candidates, and therefore assign a 
systematic uncertainty in case these form a peaking background. 

• S^! and are studied separately; the ratios are found to be consistent. We assign the 
difference as the systematic uncertainty due to detector charge asymmetry. 

• We assign an uncertainty of 1% on the efficiency for each track required to be identified as a 
kaon. 

• The uncertainty in the branching fraction of f2~ , (67.8 it 0.7)% 'W, is included. 
We make the following additional checks: 

• The desired ratio of the branching fractions is calculated from off-peak data only, and is 
measured to be 0.259 ± 0.044, consistent with the main result. 

• The data are divided up into three p* ranges: (1.8-2.7) GeV/c, (2.7-3.6) GeV/c, and (3.6- 
4.5) GeV/c; the yields and ratios calculated for each range are 0.269 ± 0.030, 0.295 ± 0.019, 
and 0.263 ± 0.031, respectively. These are consistent with being independent of p* within 
statistical uncertainties. 



5 PHYSICS RESULTS AND SUMMARY 

In summary, we observe production from the cc continuum and from T{4S) BB + X 

using the BABAR detector at SLAC. This represents the first observation of —>■ Q^K^ in B 
decays. We present a preliminary measurement of the ratio of branching fractions of to 
and S~'K~^ , determined using the cc continuum data: 

— Vt^^ — = 0.296 ± 0.018 (stat. ± 0.030 sys. . 

This represents a significant improvement on the existing value of (0.50 =b 0.21 it 0.05) 
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Table 2: Systematic uncertainties. 



Source 


Uncertainty 


Fits to mass spectrum 


0.019 


Efficiency 


0.015 


cos 9* Distribution 


0.016 


Limited Monte Carlo Statistics 


0.004 


Multiple candidates 


0.004 


Charge asymmetry 


0.001 


Particle ID 


0.006 


n~ branching fraction 


0.003 


Total systematic uncertainty 


0.030 
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